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I 

This article reviews recent theoretical and experimental 
work on a new class of topological material - topological 
Kondo insulators, which develop through the interplay of 
strong correlations and spin-orbit interactions. The his¬ 
tory of Kondo insulators is reviewed along with the theo¬ 
retical models used to describe these heavy fermion com¬ 
pounds. The Fu-Kane method of topological classification 
of insulators is used to show that hybridization between 
the conduction electrons and localized f-electrons in these 
systems gives rise to interaction-induced topological insu¬ 
lating behavior. Finally, some recent experimental results 
are discussed, which appear to confirm the theoretical pre¬ 
diction of the topological insulating behavior in Samarium 
hexaboride, where the long-standing puzzle of the residual 
low-temperature conductivity has been shown to originate 
from robust surface states. 
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1. INTRODUCTION 

Topological Kondo insulators (TKIs) are a class of narrow gap insulator in which the gap 
is created by electron correlations, but which are at the same time, topologically ordered. 
The first Kondo insulator (KI), SmBe was discovered almost fifty years ago (|T]) and today 
there are several known examples; at room temperature, these KIs are metals containing 
a dense array of magnetic moments, yet on cooling they develop a narrow gap due the 
formation of Kondo singlets which screen the local moments (H©. Topological insulators 
(TIs), discovered in just the last decade, are a new state of matter (IBV tlSI) : in which a 
band-inversion gives bulk insulator with an energy gap that is traversed by a metallic Dirac 
surface state. 

Kondo insulators were long regarded as a kind of “renormalzed silicon”, with a gap which 
is narrowed by the strong renormalizing effects of electron interactions ©■ The arrival of 
topological insulators forced a re-evaluation of this viewpoint. The large spin orbit coupling, 
and the odd-parity of the f-states led the current authors to propose Gi) that Kondo 
insulators can become topologically ordered. The recent observation of robust Gil (HI 
conducting surface states in the oldest Kondo insulator SmBe supports one of the key 
elements of this prediction, prompting a revival of interest in Kondo insulators as a new 
route for studying the interplay of strong interactions and topological order. 

In this article we review these recent developments, particularly those surrounding 
SmBe. We begin by giving a brief review of the early history and understanding of Kondo 
insulators, followed by a review of recent developments associated with these materials. 
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1.1. Key Properties of Kondo Insulators 

There are more than a dozen known Kondo insulators including the f-electron mate¬ 
rials SmBe (HJ, SmS under pressure ifT^ . YbBi 2 (l20ll . CeFe 4 Pi 2 (1^ . CeaBuPts (l22ll . 
CeRu4Sn6 ([23ll and the d-electron material FeSi l|M EH). There are also examples of 
“failed” Kondo insulators such as CeNiSn and CeRhSb (l26l \‘27ti in which the insulating gap 
appears to close in certain directions (1281129t forming a semi-metal. Kondo insulators are 
the simplest example of heavy electron materials. At high temperatures, Kondo insulators 
are simply local moment metals, with classic Curie-Weiss magnetic susceptibilities 


xiT) 


1 (M^) 

3T + 6 


( 1 ) 


that indicate the presence of a dense lattice of local moments. However, under the influ¬ 
ence of the Kondo effect, the strength of the antiferrromagnetic interaction between local 
moments and the conduction electrons grows, ultimately leading to Kondo screening of the 
local moments to produce a paramagnetic ground state. In a simple picture, the narrow 
gap of Kondo insulators reflects the energy required to break these emergent singlets. 


1.2. Early History 

Until the 1970’s magnetic materials containing f-electrons were thought to be electronically 
inert. The discovery of the first Kondo insulator SmBe by Anthony Menth, Ernest Buehler 
and Ted Geballe in 1969(|T]) changed the perspective. SmBe is a paramagnetic metal at 
room temperature, with a Curie-Weiss susceptibility characteristic of magnetic Sm^^ ions, 
yet on cooling, it evolves into a paramgnetic insulator with a tiny lOmeV gap. The discovery 
of similar behavior in pressurized SmS led Brian Maple and Dieter Wohlleben ifinilHCTll to 
propose that coherent valence fluctuations in rare-earth ions destabilize magnetism, allowing 
the f-electrons to delocalize into the conduction sea. 

Building on these ideas, Neville Mott, Chandra Varma and Yako Yafet proposed the 
idea that Kondo insulators involve a kind of excitonic ordering between localized f-electrons 
and delocalized d- electrons (I31H33II . giving rise to a hybridized band-structure with a 
gap. In another development, Sebastian Doniach introduced the the concept of a “Kondo 
lattice” dHl): a lattice of localized moments immersed in a sea of mobile electrons, described 
by the model 

^ E E ■ (2) 

i,a/3 


Here aj = {c^jpSpaC-ja) is the spin density at site j and J is the antiferromagnetic Kondo 
coupling, and the spin S = 1/2 local moment, Sj at each site j. Doniach pointed out in 
the lattice, the physics is determined by a competition between the Kondo effect, which 
tends to screen the local moments, forming Kondo singlets below the characteristic Kondo 
temperature 


Tk = Dy/jpexp 
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( 3 ) 
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and the magnetic RKKY (Rudderman Kittel Kasnya Yosida) interaction between the local 
moments, which leads to magnetic order at the characteristic temperatnre Trkky ~ J^P- 
Doniach proposed that provided Tk exceeds Trkky, the Kondo effect will overcome mag¬ 
netism, to produce a singlet gronnd-state. 

1.3. Strong Coupling model for the Kondo insulator 

A simple pictnre of the Kondo insulator is obtained by considering the Kondo lattice at 
strong coupling. Since the Kondo effect causes the coupling constant J to renormalize 
to strong coupling, the essence of the Kondo lattice can be understood by examining the 
strong coupling limit in which J is much larger than the hopping t. In this limit, the 
intersite hopping shown in 0 is a perturbation to the onsite Kondo insteraction, 

H 0{t), (4) 


and the corresponding ground-state corresponds to the formation of a spin singlet at each 
site, denoted by the wavefunction 

1^^) " n (5) 

where the double and single arrows denote the localized moment and conduction electron 
respectively, as illustrated in Fig. [^(a). 

Each singlet has a ground-state energy E = — | J per site and a singlet-triplet spin gap 
of magnitude AE = 2 J. If we remove an electron from site i, we break a Kondo singlet and 
create an unpaired spin with excited energy | J, 

Iqp"^,* t) =11^ n ^ , (6) 

while if we add an electron, the Kondo singlet is broken to create an electron quasiparticle, 
a composite involving an unpaired local moment and a doubly occupied conduction electron 
orbital 

I<1P 1* t) =lTi 0 ~ ! (7) 

as illustrated in Fig Sb)- In this fashion, the strong coupling Kondo lattice with one 
electron per site forms an insulator with a charge gap of size 3J and a spin gap of size 2J. 

If we subsequently reintroduce the hopping —t between sites, then the quasiparticles 
become mobile, as illustrated in Fig. Thus if we hole- or electron- dope the Kondo 
insulator by removing electrons, we end up with a narrow band of heavy quasiparticles or 
heavy electrons. In this way, Kondo insulators can be regarded as the parent compounds 
of heavy fermion metals. 
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Figure 1 
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Showing (a) electron and (b) hole doping of strong coupling Kondo insulator, (c) Dispersion of 
strong coupling Kondo insulator showing the formation of a heavy Fermi surface when the Kondo 
insulator is hole-doped. 


1.4. Adiabatic picture 

An alternative way to understand Kondo insulators in terms of adiabaticity was proposed 
by Richard Martin and Jim Allen l|35ll. The Kondo lattice Hamiltonian is a low energy 
limit of the Anderson lattice model, 


Halm = He + Hf + Hhyb (8) 

where 

= ^ (cLc,v+H.c) (9) 

describes the condnetion electrons, 

HatomU) = [EfUfij) + Unf^{j)nfi{j)] . ( 10 ) 

j 

describes the atomic Hamiltonian of a localized f-state at site j in energy level Ef and 
corresponding onsite Coulomb repulsive interaction U, while 

= ( 11 ) 

3 

describes the hybridization between the localized / state and the conduction electrons. At 
f/ = 0 and half filling, this model describes a simple hybridized band-structure with a 
direct hybridization gap V and an indirect gap ~ /D, where D is the half band¬ 

width, as first noted by Mott (I31II . (See Fig. (a)) By appealing to adiabaticity, Allen 
and Martin argued that as U is increased, provided neutrality is maintained, the gap will 
simply renormalize downwards. At large U, onsite charge fluctuations of the f-state can be 
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eliminated via a canonical Schrieffer WoltF transformation l|55|) and in this limit, the model 
reduces to the Kondo lattice model with J ~ /U. Adiabaticity enables one to understand 

the Kondo insulator as simply the large U cousin of the original hybridized band insulator, 
with a duality between weak and strong coupling in the two models: 

Large U ^ Small d ~ ^ 

Small U ^ Large J ~ ^ (12) 

The upshot of this discussion, is that the low energy physics of the Kondo insulator can 
be equivalently described by a renormalized Anderson lattice model, with renormalized 
parameters, V*, Ef and Ag determined by the Kondo temperature Tk (See Fig. I (b)). 
Today there are various methods for calculating these renormalizations, including path 
integral slave boson Gutzwiller and dynamical mean-field theory (ESl- 

1491) formulations of the Kondo and inhnite U Anderson models. 


(a) (b) 




U 


Figure 2 

Showing (a) hybridization of localized f and mobile d-band gives rise to a hybridization gap 
insulator with direct gap 2V and indirect gap Ag ~ jD. (b) When the interaction is turned on 
adiabatically, the band-gap renormalizes down towards the Kondo temperature. 


2. RISE OF TOPOLOGY 

The concept of topological order has its roots in the pioneering work of Robert Laughlin, 
of David Thouless, Mahito Kohmoto, Peter Nightingale, Marcel den Nijs and of Duncan 
Haldane, on the integer quantum Hall effect or Quantum Hall Insulator (I50H52II . From this 
work, the quantization of the Hall effect could be interpreted understood as a consequence 
of a topologically ordered ground-state wavefunction, in which the quantization of the Hall 
constant results from the integer-valued Chern number of the topology. Later work of 
Shuichi Murakami, Naoto Nagaosa and Shoucheng Zhang, of Charles Kane and Eugene 
Mele, and of Andrei Bernevig and Taylor Hughes, (1H1I1E31) led to the concept of the Spin 
Hall insulator: a two-dimensional topological insulator which is basically two time-reversed 
copies of the Quantum Hall Insulator. 
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In a series of remarkable papers, beginning around 2006, Joel Moore and Leon Ba- 
lents, (0, Liang Fu, Charles Kane and Eugene Mele, dn, and Rahul Roy da made the 
discovery that three-dimensional insulators can acquire a topological order through a spin- 
orbit driven band-inversion, leading to a new kind of insulator. From this new perspective, 
the presence of topological order is determined by a single Ising or Z 2 index which is posi¬ 
tive in conventional insulators (Z 2 =-|-l), but reverses sign in topological insulators (Z 2 =-l). 
Conventional insulators can be loosely considered as a miniature version of the physical 
vacuum. By contrast, topological insulators contain a internal twist to their wavefunction 
that prevents them from being adiabatically transformed into the vacuum of empty space, 
and this is why they innevitably develop surface states. We can imagine smearing out the 
surface of an insulator so that the the path from the insulator to the vacuum of empty space 
becomes an adiabatic deformation of the Hamiltonian. But if the vacua are topologically 
distinct, then the gap can not remain open along this path, or it would be possible to 
adiabatically deform the one into the other: as a result the insulating gap must collapse at 
the interface between two topologically distinct gapped states, to produce a gapless surface 
states. These states turn out to be Dirac surface states, with excitations whose spin is 
locked perpendicular to their momentum. 

Conventional surface states are incredibly volatile and rarely survive as macroscopic 
conducting surfaces, since they are highly sensitive to disorder, which leads to Anderson 
localization and surface reconstruction, which eliminates the conducting states altogether. 
However, topologically protected surface states are robust against both Anderson local¬ 
ization and surface reconstruction, leading to the unusual situation where even cracks in 
the sample are conducting! The first two dimensional topological insulators were pre¬ 
dicted in mercury cadmium-telluride (HgTe/CdTe) quantum wells by Andrei Bernevig, 
Taylor Hughes and Shou-Cheng Zhang in 2006 and were discovered in 2007 by Lau¬ 
ren Molencamp and collaborators (II3|. The first experimentally realized three dimensional 
topological insulator was bismuth antimonide (Bii_i,Sbi,) ()14ll . discovered by the groups of 
Robert Cava and Zahid Hasan in 2008. 

In 2007, Liang Fu and Charles Kane showed that if an insulator has both time reversal 
and inversion symmetry (H), the Z2 index is uniquely determined by the the parities Si„ 
of the occupied Bloch states at the high symmetry points Fi of the valence band 


Z2 = Y[^{T.) 

Ci 


-1-1 conventional insulator 
— 1 topological insulator 


(13) 


where S{Ti) = Sin is the the product of the parities of the occupied states at the high- 
symmetry points in the Brillouin zone. This beautiful index formula allows one to compute 
whether an insulator state is topological, merely by checking whether Z 2 is negative. 

Each time a band-crossing between an odd and even parity state occurs, the Z 2 index 
changes sign (see Fig. [^, and provided Z 2 = —1, there will be protected Dirac surface 
states. For a cubic insulator, with high symmetry points at the F, X, M and R points, 
this formula reduces to Z 2 ~ SfSxSmSr. (The X and M point occur three times, but 
^x,M = 5x,m). For example, if Jx = — 1 is the only negative parity, this means there have 
been three band crossings at the X points, and in general this will give rise to up to three 
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Dirac cone surface states. 



Figure 3 

Showing (a) topologically trivial band insulator with Z 2 = +1 (b) band-crossing of even and odd 
parity states at an odd number of high symmetry points leads to a topological insulator with 
Z 2 = —1- Each band crossing generates a Dirac cone of spin-momentum locked surface states. 


2 . 1 . Topology meets strong correlation 

In 2010, Maxim Dzero, Kai Sun, Victor Galitski and Piers Coleman Gi) proposed that 
Kondo insulators can form strongly interacting versions of the Z 2 topological insulator. 
The key points motivating this idea were that: 

• The spin orbit coupling of f-electrons in a Kondo insulator, of the order of O.SeV, is 
much larger than the characteristic lOmeV gap of a Kondo insulator, making these 
essentially infinite spin orbit coupled systems, ideal candidates for spin-orbit driven 
topological order. 

• f-states are odd-parity, whereas the predominantly d-band conduction bands that 
hybridize with them are even parity, so that each time there is a band-crossing between 
the two, the Z 2 index changes sign, leading to a topological insulator. 

The TKI proposal provides an appealing potential resolution of a long-standing mystery 
in the Kondo insulator SmBe, which for more than thirty years, had been known to exhibit 
a low temperature resistivity plateau (|54l I55p (see Fig.[^, which could be naturally under¬ 
stood as a consequence of topologically protected surface states II 161 156|l . In 2012, teams at 
the University of Michigan GUI and the University of California, Irvine GHD, confirmed the 
existence of robust surface states in SmBg. Most recently 2014 Xu et al. have detected 
the spin-polarized structure of the surface states in these materials that tentatively confirm 
their topological character (see discussion in section 4.2). 
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3. THEORETICAL MODEL OF TOPOLOGICAL KOISIDO INSULATORS 

In this Section we review the models, which have been recently proposed to describe various 
topologically nontrivial electronic states in Kondo insulators. We begin by introducing a 
general model for the topological Kondo insulators and discussing the conjecture of adiabatic 
connectivity between band insulators and Kondo insulators. This discussion is then followed 
by the review of the tight-binding models specific to Samarium hexaboride. 

3.1. General discussion 

The Anderson lattice model (ALM), consisting of a set of localized /-electron states, hy¬ 
bridized with conduction electrons, provides a basic description of the physics of /-electron 
materials. The ALM Hamiltonian has the form 


Halm — Hf He -i- Hhyb- 


(14) 


Here Hf describes the /-electron system and can be written as follows: 


(15) 


(f) 

where are the hopping amplitudes for the /-electrons between neighboring sites on the 
lattice, fia is a fermionic annihilation operator and [// > 0 is the strength of the Coulomb 
repulsion between the /-electrons on the same site. The interactions between the conduction 
electrons are ignored. 

In the simplest models of Kondo insulators the /-electrons are considered as dispersion- 
less, localized states, tlj = efSij. However, in the context of topological Kondo insulators, 
where the hybridization contains nodes, a more general model with nonzero off-diagonal 
elements of is required (see our discussion below). Inside a crystal, the intrinsic de¬ 
generacy of the Z = 3 /-states, = 2(2Z -I- 1) = 14 is lifted by the spin-orbit coupling and 


the crystalline fields, and the subscript a in (141 labels the components of the low-lying 
spin-orbit split /-multiplet. 

The second term in (1141) accounts for a band of non-interacting conduction electrons: 


H = fit c b 


(16) 


ij(T a,b 


where cr labels the projection of conduction electron’s spin, (a, 6) accounts for the 

degeneracy of the corresponding conduction orbitals due to non-zero angular momentum 
and is the hopping amplitude. Note that the intersite hopping terms in (15) and in 

(|16[) are not necessarily limited to the nearest neighbor sites. Lastly, the hybridization term 


in the Hamiltonian (14) is 


Hhyb — 






fLcj. 


ia,ja(T J 


(17) 
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with Viacr,ja being the hybridization matrix diagonal in momentum space: 

Viaa,j a — ^ ^ hat7,a(k) exp[ ik ■ )]. (1^) 

k 

The important point, is that all the spin-orbit coupling, all the topology is hidden inside the 
hybridization matrix. The specific form of the hybridization matrix Kio-,Qf(k) will depend 
on the multiplet structure of the f-electrons as well as on the angular momentum of the 
conduction electrons. However, typically the conduction electrons are composed of even 
parity 1 = 2 d-states, hybridizing with the 1 = 3 /-electrons, which have odd parity ®. 
It follows that at the high symmetry points in the Brillouin zone k = kr, parity is a good 
quantum number and as a result, the odd-parity f-states and the even parity d- states can 
not mix, so that H(kr) = 0. More generally, since the Hamiltonian must be even under the 
parity operation, any hybridization matrix that mixes even and odd parity electron states 
must itself be an odd-parity function of momentum 

I4a,Q(-k) = -I4a,Q(k), (19) 


i.e. the elements of the hybridization matrix are odd functions of momentum. We can 
see this more detail by noting that at a high symmetry point, the effect of reversing kr 
is the same as shifting it by a reciprocal lattice vector G, —kr = kr -I- G, and since the 
hybridization is periodic in G, i.e V (kr + G), it follows that l/(—kr) = H(kr-I-G) = '(/(kr). 
Comparing this with (19 ), it follows that V (kr) = 0 as expected. The presence of nodes in 
the hybridization at the high symmetry points kr is central to a description of topological 
Kondo insulators. 


3.1.1. Non-interacting limit and adiabatic continuity. Let us consider the simplest case 
when the interaction between the /-electrons is zero, Uf = 0. In this case, the model 
Hamiltonian ( |14[ ) describes two hybridized bands of non-interacting electrons as illustrated 
in Figure]^ (a), with a direct gap of order V and an indirect gap of the order ot\V\^/D (where 
\V\ is the typical size of the hybridization and D is the conduction electron bandwidth). 
The insulating ground state in this model corresponds to the case of half-filling, so that the 
lower bands become fully occupied. For example, in the simplest case of a single-Kramers- 
degenerate conduction electron and /-orbital Kramers doublet, an insulating state is realized 
when there is exactly one conduction electron and one /-electron per site. 

Does the system remain insulating if we adiabatically switch on the interaction [//? In 
their pioneering work Martin and Allen (p5l[58l[59)l argued using Luttinger’s theorem, that 
as long as the band remained half-filled, there would be no reason for the gap to close (See 
Fig. I^a). Numerical studies within the self-consistent fluctuation exchange approximation 
have largely confirmed this conjecture at least for moderate values of Uf. (EHml Theoretical 
schemes based on the variational wave-functions (EH also seem to suggest the validity of 
the Martin-Allen conjecture even for t// —>■ oo. 

However, topology complicates adiabaticity arguments (1161 1561) : : in particular, as we 
saw in the last section, when the parities of the conduction and f-bands are opposite, the 
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Figure 4 

(a) If we ignore the effects of topology in a conventional Kondo insulator, the interaction can be 
turned on adiabatically. When the interactions are turned on, the lower band is pushed into the 
upper band. Two bands of the same parity will always repel one-another and will not cross when 
the interactions are turned on. (b) When interactions are turned on in a topological insulator, 
they can lead to band-crossing and a topological phase transition. Here, interactions cause an 
f-band to push up into a d-band. Since the two bands have opposite parity, they do not hybridize 
at the high symmetry point so band-crossing occurs, leading to a topological phase transition. 


hybridization vanishes high symmetry points, and this opens up the possibility that interac¬ 
tions will induce band-crossing, changing the topology of the ground-state. For example, let 
us assume that in the non-interacting limit, the ALM Hamiltonian is topologically trivial, 
with a completely filled band of f-states (See Fig. [^b). For a system with time-reversal 
and inversion symmetry, the Z 2 topological invariant = 0,1 is determined by the parity 


operator eigenvalues m- For the ALM and taking into account (|19[) it is simply given by 


^2 = (-1)'' = H sign[ec(k„) - £/(k„)], 


( 20 ) 


where km is a momentum at one of the eight high-symmetry points of the 3D Brillouin zone 
and ec,/(k) is the dispersion of the conduction and the /-electrons correspondingly. As we 
switch the interaction Uf adiabatically the conduction d-band and f-bands will renormalize, 
with the f-level moving upwards relative to the conduction bands due to their stronger 
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Coulomb interaction, ec,/(k) —t £c,/(k). If as a result of this process, the f- and d- bands 
cross at an odd number of points, then a topological phase transition will take place into 
a topological insulator (see Fig. |^b). Since two topologically distinct states cannot be 
adiabatically connected without closing the gap, this phase transition must be accompanied 
by a momentary closing of the gap at a particular critical values of ?7/. In this way, 
adiabaticity can break-down. Recent studies by Werner and Assaad m of the periodic 
ALM on the two-dimensional square lattice have found that increasing the strength of the 
Hubbard interaction leads to the series of transitions between normal insulator (small Uf) 
into a strong topological insulator with the band inversion at the F point and then into 
another strong topological insulator, in which bulk bands invert at the X point of the 
two-dimensional Brillouin zone. 

From this discussion, it becomes clear that in order that the Alien-Martin idea of adia¬ 
baticity to topological Kondo insulators, then the interaction and the strength of the spin- 
orbit interaction in the Hamiltonian must be tuned so that no band crossing takes place. 
Since interactions have the effect of narrowing the bands, in practice, this will mean start¬ 
ing with a Hamiltonian with Hamiltonian with a greatly enhanced value of the spin-orbit 
coupling. 


3.1.2. Model for topological Kondo insulators with Kramers doublets. The simplest model 
for a topological Kondo insulator contains a single d-band of conduction electrons hybridized 
with a single Kramers f-doublet, with nearest neigbor hopping and hybridization. Fig. 
The resulting non-interacting Hamiltonian has the form: 


f ec(k) yd(k) 


where 


"" ^ 5 l"<5' -' v(k )") 


ic f 

ec,/(k) = —-^(cosfcj, -I- cosky + cosk^) + 


describe the dispersion of the conduction and f-band, while 

d(k) = (sin fca;, sin fcy, sin 


( 21 ) 

( 22 ) 

(23) 


is a vector in spin space, which approximates to d(k) ~ k in the vicinty of the F point. 
Notice how all the spin orbit coupling effects are held by the hybridization. In the special 
case where ec(k) = —£/(k), this Hamiltonian is a simple lattice generalization of the mean- 
field Hamiltonian for the gapped topological superfiuid phase of Helium, He-3B (I62I I63|l . 
The non-interacting band dispersion described by this Hamiltonian is simply 


^ ec(k) -fe/(k) ^ 


ec(k) - e/(k) \ 


-f F2|d(k)|2. 


(24) 


This model has been employed to great success in a number of recent publications iPP7|l. 
The low energy properties of the interacting model with the Hamiltonian 


ioil3 


(25) 
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H = (|k(T)V;a(k)(a|+h.c.) |a) = |±) 

Figure 5 

Diagrammatic representation of the hybridization process between the conduction electrons and 
/-electron Kramers doublet. 


can be analyzed by either employing the conjecture that the effect of the local correlations 
between the /-electrons leads to the renormalization of the hybridization amplitude and the 
shift of the /-energy level (see e.g. II16II56II and references therein). Alternatively, one can 
set 17/ —>■ 00 and project out the doubly occupied /-states by introducing the slave-boson 
operators ll5l 1681) . The slave-boson mean-field theory corresponds to replacing operators 
by c-numbers which are computed self-consistently. The result of the above mentioned 


procedure is a renormalization of the unperturbed bands (211: 


V ^ V = - ufV, tf ^ tf = (1 - nf)tf, 


(26) 


where n/ is an average /-level occupation number. Lastly, the position of the bare /-electron 
level €/ shifts from below the chemical potential for the conduction electrons to energies 


€f + X above it. Subsequent analysis of the topological invariant (201 yields the odd number 
of band inversions for intermediate values of Uf - mixed-valence regime - and even number 
of band inversions for n/ ~ 1 - local moment regime. (I16II56II The latter corresponds to the 
so called “weak topological insulator” with the topological invariant given by 




n 


sign[ec(k„) - £/(k„)], {a = x,y,z), 


(27) 


where Pa denotes a plane perpendicular to one of the main crystalline axes. 

An appearance of the weak topological insulating state for n/ ~ 1 raises a general ques¬ 
tion whether it would be in principle possible to stabilize the strong topological insulating 
state {v — 1) in the local moment regime. To address this issue Dzero II68|I has generalized 
conduction and /-electron Kramers doublets to N components and used the symplectic 
SP(A’) (N = 2k,k = 1,2,...) representation for the electronic operators to properly ac¬ 
count for time-reversal symmetry of the conduction and /-electron states. In agreement 
with the previous results (11611561169II he found that for N = 2 and Af = 4 there appears 
two (weak and strong) topologically non-trivial states depending on the relative position 
between the renormalized /-level and the chemical potential of the conduction band. In- 
terestengly, for the large value of > 4 there is only strong topological insulating state. 
Thus, the higher degeneracy of the /-electron multiplet favors the strong topological state. 
This result has been later confirmed by considering the tight-binding model for Samarium 
hexaboride (SmBe) with the fourfold degenerate /-orbital multiplet. (I70II 
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3.2. Tight-binding models for Samarium hexaboride 

Samarium hexaboride has recently emerged as a prototypical candidate for the first exper¬ 
imental realization of strongly correlated topological insulator. Perhaps one of the most 
intriguing problems concerns the role the electronic correlations play in determining the 
parameters of the helical surface states such as the effective mass of the surface electrons 
and characteristic length scale on which surface states penetrate into the bulk. Thus, the 
formulation of the realistic and yet tractable tight-binding model presents a first important 
step towards a better understanding of the physical properties of SmBe. Below we review 
the recent theoretical models which have been put forward to describe the formation of the 
topological surface states in this material. 

The key insight for building realistic tight-binding model for SmBe came from a detailed 
first-principle calculations of Yanase and Horima (1711) and of Antonov et al. II72II Their 
calculations show that the Samarium 4/ orbitals hybridize exclusively with Samarium 5d- 
orbitals dzll ED). The crystal field splitting of the d-multiplets leads to a low-lying Cg 
doublet at the T point. Away from the T point, the Cg orbitals split into two Kramers 
doublets, the lower one dipping down at the X point, where it dives through the the 
4/ bands. Hybridization between the two bands forces 4f states from the valence to the 
conduction band, forming heavy 4f electron band pockets at the X points. Once the d-band 
crosses through the f-band at the three X points, so long as there are no other crossings, 
the resulting non-interacting band-structure is innevitably topological, independently of the 
details of the f-nmltiplets (See Fig. [^. 


Figure 6 



Schematic illustration of the band-crossing between d- and f- states at the X point in SmBg. (a) 
Bands uncrossed. The filled 4f® band of f-electrons is a conventional insulator, (b) Bands crossed: 
the d-band cuts beneath the f-band at the A-point, displacing an odd parity f-state from the 
valence band to the conduction band. The resulting (—1)^ sign reversal in the Z 2 index gives rise 
to a topological insulator. 
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Infact, in the cubic environment, the six J = 5/2 4f orbitals of the Samarium split into 
a Ft doublet and a Fg quartet. LDA studies iniizii) suggest that the physics of the 4/ 
orbitals is governed by valence fluctuations involving electrons of the Fg quartet and the 
conduction Cg states, e~ + 4/®(Fg“^) ^ 4/®. The Fg“^ (a = 1,2) quartet consists of the 


T|): 


Fg ') = |±|). This then 


following combination of orbitals: IFg'^^) = y | |±|) + • 
leads to a simple physical picture in which the Fg quartet of /-states hybridizes with an eg 
quartet of d-states to form a Kondo insulator. 

In 2011, Takimoto II73II introduced a tight-binding model for SmBg in which the hopping 
amplitudes in the Hamiltonian (141 are non-zero for nearest- and next-nearest-neighbors, 
while hybridization involves the nearest-neighbor overlap integrals only. The values of the 
hopping amplitudes were adjusted to fit the LDA band structure results, while the effect of 
interactions between the /-electrons is modelled as a renormalization of the bare /-energy 
level and the hybridization. In Takimoto’s model, a singlet d-like orbital inverts with an 
/-like orbital at the X point of the bulk Brillouin zone, while the remaining two bands 
remain inert. This band inversion at the X points implies the existence of three Dirac 
cones on the surface: one at the surface F point and two at the X points. Interestingly, 
the corresponding Fermi velocities for the electrons at the F point are the same, while the 
Fermi velocities at the X are strongly anisotropic. m 

Alexandrov and collaborators (ESI have considered a simpler tight-binding model only 
nearest neighbor hopping amplitudes only, showing that the LDA band structure results are 
recovered by using the appropriate Slater-Koster (l74| ratio of the overlap integrals between 
d^ 2 _y 2 and d^z‘ 2_,^2 orbitals. By setting [//—>■ oo and employing the slave-boson mean-field 
theory, it was shown that strong topological insulator state extends all to way to the local 
moment regime, n/ ~ 1. It was also emphasized that quite generally, cubic topological 
insulator can only be realized when bands invert at the X or M points, since at the F and 
R points the d- and /-bands are doubly degenerate and therefore remain topologically inert 
- the parity eigenvalues at these points are always positive. 

Legner and collaborators have considered an even simpler model for cubic topological 
insulators. lIMIl . involving one Kramers degenerate d- and /-orbitals. These authors have 
focused on the general topological properties which are possible in the presence of the cubic 
symmetry. The hopping terms in the tight-binding model If64| involve amplitudes between 
nearest-, next-nearest- and next-next-nearest neighbors. As a result, for various choices 
of the hopping amplitudes various topologically non-trivial states including the one where 
d- and /-bands invert at the X points. Legner et al. have also derived the surface state 
dispersion and found that the Fermi velocity for the electrons is 


VF = 4;|I7| 


\tftd 


{td - tfY 


(28) 


This result implies that the effective mass of the surface electrons m* = pf/vf is quite 
heavy since the hybridization amplitude is small compared to other relevant energy scale, 
while the expression under the square root is of the order 0(1). 
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4. RECENT EXPERIMENTAL PROGRESS ON SmBe 

Kondo insulators, in particular SmBe are prototypical strongly-correlated mateirals featur¬ 
ing f-electron physics and some highly unusual electronic properties. As discussed earlier in 
this review, the theory of Kondo insulators has predicted in some KI materials the existence 
of a topologically nontrivial surface state with an odd number of Dirac surface bands. In 
the past two years tremendous experimental progress has been made on SmBg seeking to 
identify the existence of the metallic surface state, to test the topological properties of the 
surface state and to search for correlation physics. Below we review some of these experi¬ 
mental works. We note that despite the rapid progress of the field, questions remain about 
the exact nature of the surface state, the consequence of multiple surface bands, and un¬ 
usual surface properties due to interaction within the surface state and between the surface 
and bulk Kondo lattice. The surface state of SmBs (and possibly other TKI candidates) a 
good experimental platform to explore unusual quantum phases involving strong electron 
correlations, heavy fermion physics and topological order. 

4.1. Transport evidence for the metallic surface state 

SmBe (Fig. is a classic Kondo insulator ®. As discussed above, the strong inter¬ 
action within a localized periodic dense array of localized f-magnetic-moments leads to a 
re-organization of the electronic structure at low temperatures m resulting in an energy 
gap driven by hybridization between conduction electrons and the highly renormalized f- 
electrons. When the Fermi level lies within the Kondo energy gap, an insulator is formed. 
As a result, SmBe behaves as a correlated metal at high temperatures o, and becomes 
insulating below 40 K with the opening of Kondo gap ®, yielding an activated diverging 
resistance at low temperatures. This is illustrated in Fig. [^. However, at even lower tem¬ 
peratures below 4 K, a peculiar resistance saturation was found in the original experiments 
of Menth et al. o and Allen et al. dH (see Fig. &)• This additional metallic con¬ 
duction channel was initially suspected to originate from bulk impurity states o, but the 
improvement of the sample quality does not seem to reduce the low temperature saturation 
behavior (II1[TS1E1[751[7H1). The theory of topological Kondo insulator by Dzero et al. 
predicted the existence of a topologically protected metallic surface state (TSS) within this 
Kondo gap, naturally accounting for the resistance saturation. Experiments soon followed 
to verify whether the low temperature conduction does indeed occur on the surface. 

In 2012, signs of a low-temperature surface state were obtained while investigating the 
capacitive self-oscillation effect in SmBe single crystals. Recording the Lissajous plots of 
current-induced voltage oscillations, Kim et al., m detected an anomalous capacitance 
component (/iF) in mm-sized SmBg crystals below 4K. This was almost exactly the onset 
temperature of the saturation resistance, hinting at a common origin. The response of the 
SmBe crystal closely resembles that of a 2 H resistor in parallel with a 2/rF capacitor, as 
shown in Fig. Since a RC circuit is equivalent electronically to an insulating material 
encapsulated with a conducting surface in a bulk material, the surface conduction picture 
was a natural explanation for the anomalous capacitance. However, despite this strong 
evidence for surface conduction, alternative scenarios could not be ruled out (e.g. having 
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Figure 7 

(a), Crystal structure of SmBe . (b), Adapted from Q, activated resistance of SmBe 6 below 50 
K due to the opening of Kondo gap, (c), Adapted from llsill . resistance saturation below 4 K. 


many insulating grains embedded in a conducting bulk could also lead to an equivalent RC 
circuit, as depicted in the inset in Fig. [8f). This motivated further transport studies. 

One way to prove that SmBg is truly insulating in its bulk is to show that the bulk 
density of states (DOS) remains unchanged below the resistance saturation temperature. 
This can be done by measuring the resistance of a very fine contact to the sample: a tech¬ 
nique called Point-Contact Spectroslcopy (PCS). If the low temperature residual conduction 
arises from bulk impurity state or the above-hypothesized metal-insulator mixture instead 
of a metallic surface, the bulk electronic state would necessarily change when the coherent 
conduction channel develops at low temperatures. Zhang et al. dZSI performed PCS on 
SmBe using Ag particle point contacts. They recorded the point-contact conductance as 
functions of both Ag-SmBe voltage difference and temperature (Fig. [8f) . The former 
conductance-volt age relation, or conductance spectra, give valuable information regarding 
the bulk electronic structure at a given temperature. Zhang et al. (EHl) observed the onset 
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Figure 8 

(a), Adapted from Ref. ((T8)|. From the Lissajous voltage-current plots SmBg is found to behave 
as a resistor connected in parallel with a 2/rF capacitor. The anomalous capacitance onsets at 
4 K. (b), Adapted from II78II . temperature dependence of Point-Contact Spectra (PCS) shows the 
appearance of a zero bias feature below 100 K, indicating an opening of the Kondo gap in the 
bulk. However, the PCS remains the same as the temperature passes 4K, showing that the bulk 
remains insulating while additional metallic conductance starts to dominate transport. 


of a zero-bias conductance dip below lOOK, signaling the opening of the Kondo gap due 
to electron correlations. Perhaps the most striking Ending of the experiment is that the 
conductance spectra remained unchanged while the temperature was lowered from 10 K to 
2K, where the sample resistance saturation occurred. This proves directly that during the 
resistance saturation, the bulk of SmBg remains insulating with an almost constant Kondo 
gap. The most natural explanation of this phenomenon is to assign the low-temperature 
metallic conduction to the surface. 

The low-temperature surface conduction could also be demonstrated somewhat more 
directly by noticing that the bulk and surface contributions to conduction vary differently as 
one changes the size and shape of the sample. Wolgast et al. GZI) designed and fabricated 
a thin-plate-shaped SmBg crystal (Fig.|^) in order to enhance the difference between bulk 
and surface contributions. The electric current leads were placed on the top and bottom 
sides of the sample, which was sandwiched between two insulating silicon pieces. For the 
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Figure 9 

(a), Adapted from Ref. (Ell, transport measurement with various current-voltage configurations 
on a thin plate SmBg crystal. Depending on bulk or surface conduction, the relative size of the 
voltages would be significantly different, (b), Adapted from 1I18I I. Hall effect measurements at 
different thicknesses in a single wedge-shaped SmBg Hall bar. The Hall voltages were found to be 
inversely proportional to thickness at high temperatures, indicating bulk conduction. Below 4K, 
the Hall voltage becomes thickness independent, signaling 2-dimensional surface conduction at low 
temperatures. 
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bulk conduction in this SmBg plate, the current can zip through the middle thin part of the 
crystal; while the surface current has to travel along the surface around the corners of the 
sample. Therefore, voltage leads, placed at various locations on the sample, could register 
these two vastly different current paths. In their experiments Gil three voltages (Fig. 
are measured as a function of temperature. After comparison between theoretical simula¬ 
tions and experimental results, it was found that surface conduction indeed dominated at 
low temperatures, while bulk contribution dominated at high temperatures (nzj. 

Another transport quantity that would differ sharply between surface and bulk conduc¬ 
tion mechanisms is the Hall effect. In a perpendicular magnetic field, the transverse or Hall 
voltage is inversely proportional to the thickness of the Hall bar in the case of the usual 
bulk conduction. However, if surface conduction dominates, the thickness of the Hall bar 
is irrelevant and the Hall voltage would be thickness-independent. Therefore the thickness- 
dependence of the Hall effect should directly indicate whether bulk or surface conduction 
dominates. A technical difficulty associated with such an experiment is the innevitable 
variability in crystal thicknesses. 

This problem was solved in an experiment by Kim and Thomas et al. d) who measured 
Hall voltages at different thicknesses, within the same crystal using an SmBg Hall bar 
shaped into a thin and long wedge (Fig. [^). In smooth, well-polished SmBg wedges, they 
found that the high temperature Hall resistance Rxy = Vxy/I was inversely proportional 
to the thickness, consistent with bulk transport m . The temperature dependence (Fig. 
1^) of Rxy follows an activated behavior with a transport activation gap of 38K, confirming 
the insulating nature of the bulk. However, at temperatures below 4 K, the Hall resistance 
Rxy becomes thickness-independent, proving the domination of surface conduction. At 
these temperatures, the Rxy is also temperature independent, consistent with the metallic 
nature of the surface conduction. Without a magnetic field, the qualitative evidence for low 
temperature surface conduction was also obtained in Ref. d by comparing the local and 
non-local voltages on thin-plate samples, as was used in non-local voltage measurements 
dZSI in quantum Hall effect systems to identify the then-debated edge state conduction. 
This measurement is conceptually equivalent to the thin-plate measurement described above 
d, both taking advantage of the vastly different paths of surface and bulk conductions. 
Several samples were used for both Hall and non-local measurements, including purposely 
scrapped samples. And the surface dominated conduction was found to hold independent 
of the crystal surface direction or surface quality, which is consistent with a topologically 
protected surface state that can only be destroyed by broken-time-reversal, e.g. magnetic, 
perturbations. 

Perhaps an even more vivid demonstration of the low temperature surface conduction 
is to show that the longitudinal resistance doesn’t depend on the thickness of thin samples. 
In such an experiment, special care needs to be taken to preserve the surface quality before 
and after the thinning process. Nevertheless positive results have been demonstrated by 
several groups now (Pn 17711751 [SCTll. 
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4.2. ARPES studies of the surface state 


The electronic structure of the surface state can be probed directly by angle-resolved pho¬ 
toemission spectroscopy (ARPES). An observation of odd number of Dirac surface bands 
would be a direct confirmation of the topological nature of the surface state. In addition, 
spin-resolved ARPES can in principle investigate the spin structure of the surface state. 
This has led to a flurry of recent ARPES studies of SmBe . Because of the limited energy 
resolution of ARPES relative to the small Kondo gap of SmBe , the detailed spin structure 
of the surface states is currently still contraversial (l57l 1811) . Important insights regarding 
the energetics of the surface state have been nevertheless been obtained by several recent 
ARPES studies. Before we review these experimental results, it would be helpful to first 
take a look at the predicted SmBe energy spectrum for a topological Kondo insulator 
obtained by Lu et al. II82I) . Another calculation has been carried out by Alexandrov et 
al. (I70II . One calculated energy spectrum II82I) is shown in Fig. 10 1 for the first Brillouin 


zone (BZ) in (001) surface, where high-symmetry points are labeled as P, X, and M. The 
blue lines are bulk energy levels with a Kondo gap at the Fermi level. And the red lines are 
surface states, showing three Direct cones at P, and the two X points. 


The ARPES results of Neupane et al. (18311 are replotted in Fig. 10 r for direct compar¬ 
ison. The topology of the measured Fermi surface at 6 K agrees well with the topological 
Kondo insulator calculation (Fig. |10^ ). As predicted, three low-lying metallic states near 
the Fermi energy exist at the P, and the two X points. The two X pockets are quite large, 
suggesting that the Fermi energy is far away from the (possible) Dirac points at X points. 
Two other important results are obtained by Neupane et al. dHSI from the temperature 


evolution of the ARPES spectral intensity (Fig. 10 r). The first is that the in-gap metallic 


state only emerges at temperatures below 15 K for both F and X pockets, suggesting the 
nontrivial origin of the metallic state. The second finding is that the metallic in-gap state 
does not disappear after a 6 to 50 K thermal recycling, which demonstrates that the in-gap 
states are robust and protected against thermal recycling. 

Two additional ARPES results by Xu et al. (18411 and Jiang et al. II85D were reported 
at the same time which are consistent with the topological Kondo insulator scenario for 
SmBe. In particular in Ref. (|84l). Xu et al. extracted the energy dispersions of the P band 
(X pocket) (Fig. |11^ ). The linear dispersions demonstrate the two-dimensional nature of 
the pocket. Unfortunately the Dirac points could not be observed clearly in this study. For 
the p band, Xu et al. found that the intensity diminishes suddenly at 20 meV below Fermi 
level, corresponding to the hybridization gap edge between f and d electrons. They argue 
that this may prohibit observing the Dirac point formed by the bands crossing each other. 
In ()85|l Jiang et al. performed circular dichorism (CD) ARPES measurements. CD ARPES 
may shed some light on the spin structure of the metallic state. Jiang et al. calculated the 
differential map of the photoemission intensities with right-circularly-polarized (RCP) and 
left-circularly-polarized (LCP) lights as repotted in Fig. The anti-symmetric pattern 
with respect to the F-X line resembles that of a topological surface state. However, it is 
not yet sufficient to make conclusive statements about the surface spin texture or Berry’s 
phase (I83II due to several complexities related to CD ARPES. 

Parallel ARPES studies on SmBe were carried out by Zhu et al. II86II and Denlinger 
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Figure 10 

(a), Adapted from Ref. calculated electronic structure of SmBg for the first Brillouin zone 

of the (001) surface (inset). The blue lines are bulk energy levels and the red lines are surface 
states, showing three Direct cones, (b), Adapted from lEg, topology of the Fermi surface 
measured by ARPES, consistent qualitatively with the theoretical calculation (I82II . The left 
figure is a temperature evolution of the ARPES spectral intensity, showing the onset of the in-gap 
(surface) state at low temperatures. 


et al. (I87II respectively. In Ref. (18611 . Zhu et al. propose an alternative explanation for 
the conducting surface state, suggesting that it originates from boron dangling bonds on 
the (001) crystal surface, i.e. that it is a non-topological polar surface. In (I87II Denlinger 
et al. found high ARPES intensity at the H points, where the momentum-location of the 
H-point along the P-M direction coincided with the polar metallic surface state claimed by 
Zhu et al. (|86|l. However, Denlinger et al. (|87|1, found that the H-point is gapped at low 
T temperature, suggesting that the polar surface state found by Zhu et al. dMI is in fact 
insulating and hence unrelated to the metallic surface state. Due to the limited resolution, 
it is somewhat hard to judge which scenario is correct at the moment. However, the (001) 
polar metallic surface picture by Zhu et al. (I86II seems to contradict available transport 
experiments that showed the existence of surface metallic conduction on both (001) and 
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Figure 11 

(a), Adapted from Ref. (I84ll . extracted dispersions of the ^ band (X pocket) for different photon 
energies. The linear dispersions demonstrate the two-dimensional nature of the pocket, (b), 
Adapted from ll85ll . the differential map of the photoemission intensities with 

left-circularly-polarized and right-circularly-polarized light. The anti-symmetric pattern resembles 
that of a topological surface state. 


(Oil) surfaces. The latter surface is not polar and therefore can not host a polar metallic 
surface state. 

Finally we mention the current status of spin-polarized ARPES measurements. In 
weakly interacting topological insulators, the classic test for topological Dirac states, is to 
directly observe their spin texture using spin-resolved ARPES measurements ll88|l . However, 
in SmBg, the narrowness of the bulk insulating gap makes this a highly challenging mea¬ 
surement, the results of which are still contraversial. The first spin-polarized ARPES mea¬ 
surements have been carried out on SmBe were carried out at the Swiss Light Source (I57II . 
These measurements show of spin-momentum locked states on the 100 surfaces in the vicin¬ 
ity of the X points in the surface Brillouin zone (see Fig. 121, supporting the presence 
of topological Dirac states. However, spin-polarized ARPES measurements carried out at 
the Helmholtz-Zentrum, Berlin (I81II have led this group to question the topological inter¬ 
pretation. In particular, their measurements around the P points appears to support the 
presence of Rashba-split, conventional surface states of a trivial insulator. The resolution 
of this experimental controversy awaits improved resolution and further measurements. 
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Figure 12 

(a) Adapted from Ref. 153, showing the spin polarization of the surface states observed along the 
high-symmetry lines. The red lines indicate the locations of spin measurements in the kx, ky 
plane, labelled as Cl, C2, C3 and C4. The arrows indicate the measured spin polarizations at the 
positions of the black dots. In (b). The low-energy excitations along the high-symmetry line 
M — X — M. (c)- (e) show the spin polarization measured at 30, 30 and 26 eV respectively at 
points Cl, C2 and C3, with the spin polarization along the x, y and x directions respectively. 


4.3. STM studies of the surface state 

Scanning Tunneling Microscopy (STM), is another useful surface probe with the potential 
to reveal the properties of the surface states with atomic spatial resolution. Spectroscopic 
measurement of the tunneling conductance as a function of bias voltage, or Scanning Tun¬ 
neling Spectroscopy (STM) further reveals the local density of states near the Fermi energy. 
Therefore a temperature-dependent STS study would tell quite a lot about the energetics 
of bulk and surface state, complementary to ARPES measurements. 

STM and STS study on the (001) surface SmBe by Yee et al. (I89II revealed that the 
surface of SmBe is quite complex and far from an ideal. Their STS morphology image of 
the normal (001) surfaces shows that much of the surface reconstructs to a 1 x 2 structure in 
which half the surface is samarium-terminated (See: Fig. |13^). Other regions of the surface 
appear amorphous and only small regions involve a pristine termination. It is unclear at this 
stage if the reconstruction is due to the particular cleavage method or whether it is intrinsic 
to SmBe . Nevertheless, this STM study suggests that the ARPES studies, which have 
optical spot sizes of hundreds of microns could have picked up a mixture of spectra from a 
variety of different surface morphologies. This may explain in part the limited resolution of 
SmBe ARPES measurements. 

In the STM study by Yee et al. (I89II a common STS spectra feature was observed 
for both 1x1 and 1x2 regions: a minimum of dl/dV at the Fermi level and a gap 
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Figure 13 

(a), Adapted from Ref. (l8^ . STM morphology image of a 2 X 1 half-Sm terminated surface and 
its STS spectra for various temperatures, (b), Adapted from Ref. (I90II . STM morphology image of 
one of the four phases found by the STM study of Duan et al. The temperature dependence of 
STS spectra for the “donut” D2 phase. 


(dl/dV peak) slightly below this minimum in energy. Fig. 13 1 shows the STS spectra for 
the case of 1 x 2 region, where the peak is centered at - 8 meV and becomes smaller at 
higher temperatures. Extrapolating to high temperature, Yee et al. argued that the peak 
(gap) would vanish at 40K, suggesting its hybridization nature. They also observed residual 
spectral weight spanning the hybridization gap down to the lowest temperatures, which is 
consistent with a topological surface state. Qualitatively similar findings were observed for 
the 1x1 surface region, although the energy is shifted by tens of eV. This testifies the 
robustness of the hybridization gap. In the other STM work ([90ll by Ruan et al., four 
types of surface morphology were found, highlighting again the complexity of the SmBe 
surface. The STM image for the “Donut” D2 phase is shown in Fig. [T^d with a cartoon 
illustrating the “Donut” structure. The temperature dependence of the STS spectra (Fig. 


13 a) is qualitatively similar to those found by Yee et al. The data revealed the emergence 
of a resonance peak below 40K, signaling the opening of the Kondo gap. Summarizing 
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Figure 14 

(a), Adapted from Ref. m, contrasting behavior between non-magnetic Y-doped and magnetic 
Gd-doped SmBg samples. Magnetic doping destroys the metallic surface state, (b), Adapted 
from Ref (H, de Haas van Alphen quantum oscillation from a SmBe sample. Extrapolation of 
the Landau level index may indicate a Berry phase of tt. (c), Adapted from Ref. ll93t . weak 
anti-localization effect in SmBg as expected from a topological surface state. 


the body of all STM works, the large zero bias conductance that persists to the lowest 
temperature may come from the topological surface state. However a smoking gun STM 
evidence for TKI is still lacking: the quasiparticle interference patterns showing surface 
state with helical spin texture. 
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4.4. Is the surface state topological? 

The experiments reviewed up to this point, provide strong evidence for the existence of 
the metallic surface state that emerges from an insulating bulk state at low temperatures. 
The metallic surface state could be conveniently probed by transport methods due to the 
insulating nature of the bulk. Although the low temperature emergence of the surface state 
(PI) and the antisymmetric ARPES pattern with opposite chirality of light (Pl> seem to be 
consistent with a topological surface state, more direct evidence is still needed to tell whether 
the surface state is really topological. In an ideal world, energy- and spin-resolved ARPES 
conducted at low temperatures below the resistance saturation would confirm whether the 
spin texture of the surface state has a Berry phase. However, the required high energy 
resolution (few meV) is beyond the current capability of both synchrotron ARPES (40 meV) 
and laser ARPES (10 meV). Other difficulties derive from the limited low temperature 
range and spin resolution with spin-resolved ARPES techniques. It is noted that other 
TKI candidates may possess larger Kondo gaps, making them to be within the resolution 
of current spin-resolved ARPES capabilities. As mentioned above a smoking gun STM 
evidence for the topologic surface state would be the quasiparticle interference patterns 
showing surface state with helical spin texture. This might be within the capability of 
current STM techniques, if a better surface without reconstruction could be prepared for 
STM studies. 

Since the bulk of SmBg is insulating at low temperatures, it is possible to test the 
topological properties of the surface state, giving partial evidence for the topological surface 
state. A topological surface state has three aspects of topological protection. First, their 
fundamental Z2 topology preserves a gapless surface state unless time reversal symmetry 
(TRS) is broken. Second, helical spin polarization prevents momentum backscattering from 
—k to k by non-magnetic impurities. Finally, the Berry phase protects the surface state 
from weak localization through time reversed paths. These collectively provide a robust 
surface state with TRS conservation. Therefore for a topological surface state, it could only 
be made insulating by broken TRS such as those created by magnetic perturbations. In 
addition, such a surface state would become less conductive when a small magnetic field 
is applied, destroying the negative interference between time-reversed carrier paths. This 
is called the weak-antilocalization (WAL) effect (I94p . These two experimental tests have 
been carried out by Kim et al. dm and Thomas et al. iP|) and both have yielded positive 
results in agreement with a topological surface state. 

Kim et al. dm made a comparison study between SmBe crystals with magnetic 
(Gd) impurity doping and non-magnetic (Y or Yb) impurity dopings. In those thin-plate 
shaped samples, surface vs bulk conduction could be distinguished by performing thickness 
dependence measurements; when surface conduction dominates, the resistance ratio should 
be irrelevant to thickness. Kim et al. found that Yb and Y doped SmBg exhibit bulk 
conduction at high temperature and change to surface conduction at low temperatures, 
similar to pure SmBe . In contrast, the resistance of Gd doped SmBe remains inversely 
proportional to sample thickness for all temperatures, indicating bulk-dominated conduction 
even at low temperature. A perhaps more direct illustration of the contrasting impact from 
magnetic and non-magnetic dopings is plotted in Fig. |13^. Below 4 K, the resistance of 
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the Y doped sample saturates just like in pure SmBg , while the Gd doped sample shows 
insulating behavior. Only the broken TRS magnetic doping destroys the metallic surface 
state. It is noted that in this study, the amount (3 percent) of Gd doping is sufficiently 
small to not introduce order magnetic states, as evidenced from the magnetic measurements 

m- 

Weak antilocalization (WAL) II94|I is expected in a topological surface state due to 
an unusual Berry phase of tt ( 00 , which causes destructive interference between time- 
reversed electron paths and lowers the sample resistance. This effect is destroyed by a time- 
reversal-symmetry-breaking magnetic field, giving rise to magneto-resistance dip around 
zero field. Unfortunately, WAL alone can not provide conclusive evidence of a topological 
surface state because strong spin-orbit coupling can also produce a WAL effect. Indeed 
Thomas et al. observed clear WAL effect in SmBe samples at low temperatures (15311 . As 
shown in Fig. |13| :, the WAL manifests itself as a conductance peak at zero field, with the 
height of the peak decreases at higher temperatures. The shape of the WAL feature fits 
well to the Hikami-Larkin-Nagaoka (HLN) equation (I95II . The temperature dependence of 
the coherence length could be extracted using the fittings, suggesting that below 80 mK, 
electron-electron interaction dominates the scattering of transport current. The WAL effect 
has also been detected by Nakajima et al. (1961) . 

The topological aspect of the surface state could also be probed by quantum oscillations: 
the periodic variation of physical properties of the sample as the Landau levels pass through 
the Fermi level in a changing magnetic field. For a topological surface state with odd number 
of Dirac points, the Berry’s phase would be tt instead of the usual value of 0. As a result, 
the Landau level index would be half-integer at the large field limit. Using a sensitive 
torque magnetometry setup, Li et al. dSH) have recently reported the detection of de Haas 
van Alphen oscillations (dHvA) in the magnetization of SmBe crystals. Fig. |13^ shows 
two types of dHvA oscillations in that sample, identified as the (001) and (Oil) surfaces 
respectively. The Landau level index of both oscillations extrapolate to -1/2 at the high 
field limit, suggesting a non-trial Berry phase of tt in both energy bands. Their experiment 
finds that the frequencies of the observed oscillations only depend on the perpendicular 
component of the field, as expected for a surface state. One of the intriguing aspects of the 
measurements, is that the quasiparticle effective mass obtained the temperature dependence 
of the signal is about O.lme. This small mass is unexpected on the basis of the “Heavy 
Fermion” nature of SmBg . More experimental and theoretical investigations are needed to 
understand this discrepancy. 


5. OPEN QUESTIONS 

Topological Kondo insulators represent a collision of three areas of research: topological 
matter, strongly correlated systems and surface science, and experts from each of these 
fields find themselves outside their comfort zone. It is against this backdrop that SmBe 
as a candidate TKI poses a immense new theoretical, experimental, and materials science 
challenges. 

The physics of Kondo insulator SmBe originates from a periodic dense lattice (Kondo 
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Figure 15 

(a), Adapted from Ref. dH, very low temperature resistance saturation of several SmBe 
samples, showing log(T) divergence before saturation, (b), Adapted from Ref. (I96II . hysteretic 
magnetoresistance. The difference of conductance between hysteresis lie within the range of 
quantum conductance. This could be explained by one-dimensional edge states one the surface 
between spin polarized domains. 


lattice) Q of localized magnetic moments, and is quite different from the Kondo effect 
([^ of isolated magnetic moments. Moreover, while it may be adiabatically connected to 
a topological band insulator, the system also lies at the brink of magnetism. Indeed, high 
pressure experiments show that at pressures of about 9GPa, the hybridization gap of SmBe 
closes, leading to the development of magnetic order ll98[ I99II . This raises many interesting 
questions and possibilities, for instance: 

• What is the effect of the surface on the Kondo effect? 

• Could the bulk and its excitations be different in any way from a conventional insu¬ 
lator? 

• Will magnetism emerge on the surface, and if so, what kinds of highly correlated and 
competing order are realized near the point of instability? 

• Finally, what is the effect of disorder, such as vacancies or “Kondo holes” at the 
Samarium sites? 

In a study by Thomas et al. (|93| who followed the saturated resistivity of SmBe to 
mK temperatures, the resistivity was found to exhibit a log{T) temperature dependence 
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(Fig. 15 i) before the final saturation below 40 mK. The size of the log{T) behavior is both 
sample-dependent and inconsistent with the WAL effect alone even with the correction of 
electron correlation. However, the temperature dependence can be perfectly explained by 
assuming the Kondo effect in the surface state driven by by Kondo holes, with a universal 
Kondo temperature of 0.6 K. In fact, the data from several samples could be scaled into 
a single numerical renormalization group calculation of Kondo effect behavior assuming a 
single Kondo temperature (Fig. |15^). 

Nakajima et al. reported hysteresis in magneto-resistance l(96ll in a SmBe sample, as 
shown in Fig. The hysteresis may originate from ferromagnetic state on the surface, 

which they hypothesize to arise due to Kondo holes. The edge state between magnetic do¬ 
mains of opposite chirality could in principle yield quantized (Hall) conductance. Nakajima 
et al. observed that the hysteretic part of the longitudinal conductance lies within the range 
of the quantum conductance (Fig. 15 a), hinting the aforementioned mechanism. 

It is clear from these examples, that the effects of interactions on topological surface 
states pose a major challenge. There have been a number of recent theoretical publications 
that have started to explore these issues II65I I10nlll02|l . Bitan Roy and collaborators (|101ll 
have emphasized an itinerant description of the interacting surface Dirac cones and in lllU2ll 
predict a proclivity towards excitonic and nematic instabilities. The papers (iBsiirnuii discuss 
the physics in terms of a surface Kondo lattice involving local moments interacting with 
chiral surface states. Alshuler and Aleiner have considered a simplified one-dimensional 
surface state, and show that interactions with the local moments will localize surface states. 
Alexandrov, Erten and Coleman (|65|l have proposed a similar two dimensional model and 
propose that the surfaces of topological Kondo insulators may develop quantum critical 
ground-states. 

One of the outstanding puzzles concerns the very high group velocities of the sur¬ 
face states of SmBg measured in experiment. Both quantum oscillation ll92ll and ARPES 
studies (I83H85II . show the surface quasiparticles are light. ARPES measurements indicate 
surface quasiparticles with Fermi velocities ranging from 220meVA (l85ll to 300 meVA (l83ll . 
Velocities obtained from dHvA measurements are two orders of magnitude larger than these 
values (I92II . By contrast, current theories 117011821 [10111103ll predict heavy Dirac quasiparti¬ 
cles with velocities vs ~ 30 — 50meVA. Alexandrov, Coleman and Erten (|65|) have recently 
suggested a link between these discrepancies and a possible break-down of the Kondo effect 
on the surface. They argue that the reduced co-ordination numbers of magnetic ions on the 
surface, the Kondo temperature is suppressed, leading to surface Kondo breakdown. The 
re-emergence of unscreened local moments at the surface is likely to favor the development 
of surface magnetism. Alexandrov et al. find that the release of d-electrons from Kondo 
singlets at the surface does not destroy the surface states, but instead has the effect of 
doping the surface Dirac cones and driving the Dirac point into the continuum; surface 
Kondo breakdown also increases the conduction band character of the topological surface 
states, leading to an approximately ten fold increase in the quasiparticle velocities, a result 
that may account for the high surface velocities seen in ARPES measurements on SmBe. 

The history of SmBe studies has invariably shown that new physics accompanies higher 
quality samples. Currently, SmBe crystals are usually grown using either the aluminum 
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fliix method or the floating zone method. The majority of the reviewed experiments have 
been performed on flux-grown-samples, while high quality floating-zone-grown samples also 
showed resistance saturation as reported by Hatnean et al. lITOill and numerous older 
reports. However, the materials aspect of SmBe is far from trivial. Phelan et al. (|80l) 
recently report that some floating-zone-grown SmBg crystals remain insulating at the lowest 
temperature measured, while additional non-magnetic carbon doping restores the surface 
metallic state. While the exact role of carbon doping is unclear at this stage, these result 
clearly highlight the complexity of surface chemistry in SmBe and the vast number of 
possibilities of controlling the topological and non-topological properties of SmBe . 
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